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^ ; Abstract 

Event-by-event hydrodynamics (or hydrodynamics with fluctuating initial conditions) has been 
. developed in the past few years. Here we discuss how it may help to understand the various 

| structures observed in two-particle correlations. 
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^ § 1. 3+1 hydrodynamics: comparison with data 

One of the most striking results in relativistic heavy ion collisions is the existence of struc- 
tures in the two-particle correlations (llfl&BElBB] plotted as function of the pseudorapidity 
difference A77 and the angular spacing A0. The so-called ridge has a narrow Acp located around 
zero and a long At/ extent. The other structure located opposite has a single or double hump in 
A0, its Atj extent is not well established. 

It has been suggested that the combined effect of longitudinal high energy density tubes 
(leftover from initial particle collisions) and transverse expansion is responsible for the ridge 
hydrodynamics based explanation is attractive because of the various simi- 
qq . larities between bulk matter and ridge (transverse momentum spectra, bary on/meson ratio, etc). 

' In order to compare with the ridge data, 3+1 hydrodynamics must be used. Traditionally, smooth 

initial conditions have been used (see e.g. 11131 Q). However some groups have started to 
develop an event-by-event approach: for each collision, some initial conditions are generated, 
hydrodynamics is run and results are stored. This is done many times thus mimicking expe- 
rience. This method was pioneered by the Brazilian collaboration SPheRIO Qil Il7h (some 
typical results can be seen in [H [III IIS Eh)- It nas been used also by the German group with H. 
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Petersen/M. Bleicher et al. [22]. More recently K. Werner and collaborators have also obtained 
results Q. 

The SPheRIO group uses the NeXus code TO2411 to generate initial conditions and the Smooth 
Particle Hydrodynamics method to solve the perfect fluid conservation equations. Results on 
two-particle correlations can be found in l^lllllllH]- H. Petersen/M. Bleicher et al. start with 



initial conditions generated with UrQMD and solve the perfect fluid hydrodynamics equations 
with a standard (grid) method. They have no results on two-particle correlations at RHIC. K. 
Werner et al. use EPOS initial conditions and also a standard (grid) method to solve the ideal 



hydrodynamics equations. Their results on two-particle correlations can be found in 112311 . 
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Figure 1 : NeXSPheRIO central Au+Au collisions at 200 GeV A: initial energy density and resulting two-particle corre- 
lation (p t rig > 2.5 GeV x p a t ssoc > 1.5 GeV). 



In the NeXSPheRIO approach, initial conditions have tubular structures and two particle 
correlations exhibit near- and away -side ridges as can be seen in Fig. [TJ As a double check, two 
different methods for elliptic flow subtraction have been used: ZYAM and event plane alignment. 
In the calculation by K. Werner et al., initial conditions also have tubular structures, a plane 
alignment method for elliptic flow subtraction is used and two-particle correlations exhibit near- 
side [23] and small away-side ridges 112711 . 

In addition to reproducing both the near and away- side structures, NeXSPheRIO leads to 
good qualitative agreement with various data. For example, 1) for fixed pf 18 and increasing 
passoc^ near- side and away- side peaks decrease as seen in figure [2] for central collisions (this 
is generally expected since the number of associated particles decreases). On the other side, for 
fixed p a t ssoc and increasing p tr t 18 , the peaks increase. As shown below (see Fig. 0, this is because 

the single-particle angular distribution 
becomes sharper, reducing the smear- 
ing effect when integrated over the trig- 
ger angle. This behavior is in agree- 
ment with data 1281. 2) When going 
from central to peripheral collisions, 
the near- side ridge decreases and the 
away- side ridge changes from double 
to single hump, as seen in figu re [3] 
This is in agreement with data 11280 . 
3) The correlation can also be stud- 
ied as a function of the trigger-particle 
angle with relation to the event plane. 
In Figure |4] for a mid-central window, 
the away-side ridge changes from sin- 
gle peak for in-plane trigger to double 
peak for out-of-plane trigger. For cen- 
tral collisions (not shown), it is always 
double-peaked. This is in agreement 
with data iH. 
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Figure 2: p t dependence of NeXSPheRIO two-particle corre- 
lations for 0-20% Au+Au collisions. Left column in legend is 
trigger momentum and right, associated particle momentum. 
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Figure 3: NeXSPheRIO two-particle correlations for central (a) to peripheral (f) collisions. 




Figure 4: NeXSPheRIO (20-30%) results for in-plane (a) to out-of-plane (f) trigger. 



2. 2+1 hydrodynamics: one tube model 

When using NeXSPheRIO, it is not clear how the various structures in the two particle cor- 
relations are generated. To investigate this, we use a simplified model. Fig. [T](left and center) 
shows a typical example of initial conditions (initial energy density) obtained in NeXus with 
various tubular structures along the collision axis. In the simplified model, only one of the high 
energy tubes from NeXus (chosen close to the border) is considered and the complex back- 
ground is smoothed out by using the average over many events (for details see B25I1 ). Transverse 
expansion is computed numerically while longitudinal expansion is assumed boost-invariant, un- 
til freeze out at some constant temperature. The resulting single particle angular distribution has 
two peaks located on both sides of the position of the tube placed at = as seen in Fig. \5\ 
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Figure 5: Angular distributions of particles in some different pj intervals, in the simplified model and resulting two- 
particle correlations. 



(left). This double peak structure is observed for all transverse momenta at more or less the same 
position. 

This two-peak emission is in contrast with what happens when a blast wave is assumed, 
namely the fact that high-energy tubes emit in a single direction. However its occurrence can 
be understood from Fig. [6] As time goes on, as a consequence of the tube expansion, a hole 
appears at the location of the high-energy tube (as in [30D). This hole is surrounded by matter 
that piles up in a roughly semi-circular cliff of high energy density matter, guiding the flow of 
the background matter into two well-defined directions. The outer part of the circle does not 
appear because the vacuum does not offer any resistance against the tube expansion and this part 
of the fluid is emitted isotropically. The two extremities of the cliff emit more fast particles than 
the background, this gives rise to the two-peaks in the single-particle angular distribution. The 
emission is not quite radial as shown by Fig. indicating that there was a deflection of the 
background flow caused by the pressure put by the high-energy tube. 

From Fig. [5] left, we can guess how the two-particle angular correlation will be. The trigger 
particle is more likely to be in one of the two peaks. We first choose the left-hand side peak. The 
associated particle is more likely to be also in this peak i.e. with A0 = or in the right-hand side 
peak with A0 ~ +2. If we choose the trigger particle in the right-hand side peak, the associated 
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Figure 6: Temporal evolution of energy density for the simplified model (times: 1.0, 2.9 and 6.6 fm). Arrows indicate 
fluid velocity on the freeze out surface, thicker curve labeled by the freeze out temperature 0.14 GeV 
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particle is more likely to be also in this peak i.e. with A(p = or in the left-hand side peak with 

Acp 2. So the final two particle angular correlation must have a large central peak at A0 = 

and two smaller peaks respectively at A0 ~ ±2. Fig. \5\ (right) shows that this is indeed the 
case. The peak at A0 = corresponds to the near- side ridge and the peaks at A0 ~ ±2 form the 
double-hump ridge. We have checked that this structure is robust by studying the effect of the 
height and shape of the background, initial velocity, height, radius and location of the tube l25ll . 

The results above were obtained for central collisions. We have also generalized this one 
tube model to non-central collisions. In Fig. H the single-particle angular distribution is shown 
with and without tube. When the tube is present, the single particle angular distribution may also 
have two peaks with separation A0 ~ 2 or this structure may be more hidden by the elliptic flow 
depending on the tube position. After subtracting the elliptic flow and averaging over the tube 
angular position, the two-particle correlation has a near-side peak and an away-side substructure 
that is indeed double-peaked for more central collisions and single-peak for more peripheral 
ones. We have also studied the dependence on the trigger angle with respect to the reaction plane 




Figure 7: Tube model for a non-central collision (20-30% centrality). From left to right: Energy density distribution, 
single-particle angular distribution for tube in the reaction plane, idem but tube at 60° . 



and found the correct behavior (single peak 
one). This result is shown in Fig. [8j 
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Figure 8: Two-particle correlation for different trigger 
(not tube) angle. 



in-plane trigger and double peak for out-of-plane 

With these information, we can discuss 
what happens in a more complex event such 
as a NeXus event. In such an event, only the 
outer tubes need to be considered. The shape 
of the two-particle correlations for a single 
tube (in particular the peak spacing) is rela- 
tively independent of its features so the var- 
ious tubes will contribute with rather similar 
two-peaks emission pattern at various angles 
in the single-particle angular distribution. 
For this single event, the two-particle correla- 
tion has a well-defined main structure similar 
to that of a single tube (Fig. 0) surrounded 
by several other peaks and depressions due 
to trigger and associated particles com- 
ing from different tubes. These additional 



peaks and depressions have positions depending on the angle of the tubes between them. When 
averaging over many randomly fluctuating events these interference terms disappear and only the 
main one-tube like structure is left. 



3. Conclusion 



To summarize, a unified picture for the structures observed in two particle correlations at low 
to moderate transverse momentum was presented. It is based on the presence of longitudinal 
high energy density tubes in the initial conditions. These tubes are leftover from the initial 
particle interactions. During the hydrodynamical evolution of the fluid, pressed by the strong 
expansion of the tubes located close to the border, matter (coming mainly from the background) 
is deflected into two symmetrical directions, leading to two-particle correlations with a near- side 
and a double hump away- side ridges. 

Various models have been suggested to explain the structures in the two particle correlations 
and even those based on hydrodynamics have strong differences among them: for each tube, there 
may be no emission ll30ll . emission in one direction 

Ml 111 H2] or two directions as advocated 
here. Three-particle correlations have been suggested as a possible test to distinguish between 
various scenarios. For example, three-particle correlations can be plotted as function of the 
angles A0i between the first associated particle and the trigger and A02 between the second 
associated particle and the trigger. The Mach cone model 113 ill or the AMPT model (13211 would 
lead to off diagonal peaks in this plot while certain models might not. Experimentally, such 
off diagonal peaks do appear in three-particle correlations 113 3n (see also the approach in B4II ). 
Naively, in our case, one might expect no off diagonal peaks because there are only two peaks 
in the single particle angular distribution (Fig. [9] left). However if one carefully consider the 
contributions coming from different 3 -dimensional integration domains (in <p t , <pi and 02), there 
appear also two off-diagonal local maxima in the three-particle correlation. This can be seen in 
the computed result in Fig. [9] right. So this type of correlation does not look promising to test 
our model. 

Another possibility would 
be 2+1 correlations, where 
the trigger and first associ- 
ated particle are chosen in 
different peaks of the sin- 
gle particle angular distribu- 
tion (A0i ~ 2) and the sec- 
ond associated particle posi- 
tion varies. When this cor- 
relation is plotted as function 
of A02 and A772 (pseudora- 
pidity difference with respect 
to the trigger), there should 
appear two main structures 
elongated in A772 correspond- 
ing to the two possibilities for 

the second associated particle: it comes from the same peak as the trigger or same peak as the 
first associated particle. Though it is not ruled out, it seems difficult that jet based models lead to 
this. 
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Figure 9: Three-particle correlation in A(p\-A(p2- naive expectation for the 
one-tube model and computed raw correlation. 



To conclude: usually, the initial conditions in the hydrodynamic description of relativistic 
nuclear collisions are assumed to be smooth. It seems however that each time more data require 
a knowledge of the event-by-event fluctuating initial conditions rather than smooth initial condi- 
tions (for example fluctuations in elliptic flow l35ll or perhaps the very behavior of elliptic flow 
as function of pseudorapity 113611 ). Now two-particle correlations might also offer us a chance to 
get a glimpse of the initial conditions and learn about the strong interaction. 
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